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(54) Radio environment analysis apparatus 

(57) A radio environment analysis apparatus 
includes a receiver and signal processor. The receiver 
receives a signal modulated by phase shift keying 
(PSK) using a pseudo noise (PN) code sequence by an 
array antenna at a plurality of fixed points on the plane, 
converts the received signals into intermediate frequen- 
cies or demodulated signals, and outputs them. The sig- 
nal processor processes a signal output from the 
receiver to analyze the arrival angle, delay time, and rel- 
ative power of the signal incoming to the receiver. The 
signal processor has a plurality of normalization units, 
first estimation unit, arithmetic unit, and second estima- 
tion unit. The normalization units individually normalize 
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signals from the receiver and output them as a normal- 
ized signal group. The first estimation unit calculates a 
covariance matrix, and estimates the arrival angle of 
each incoming signal. The arithmetic unit calculates a 
weight of the array antenna for suppressing incoming 
signals except for an incoming signal having the esti- 
mated arrival angle, and calculates the product of the 
weight of the array antenna and the normalized signal 
group. The second estimation unit calculates a covari- 
ance matrix and estimates the delay time and relative 
power. 
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Description 

[0001] The present invention relates to a radio environment analysis apparatus and. more particularly, to a radio envi- 
ronment analysis apparatus for measuring the arrival angle, delay time, and relative power of a signal incoming to the 

5 antenna in a multipath environment such as an urt^an or indoor space. 

[0002] To realize high-speed digital communication, frequency resources must be effectively used. Many studies have 
been made for techniques such as a method using the milliwave band where a wide band is essentially ensured from 
a specific band, a method of reducing a repetitively used frequency with a decrease in cell area, and dynamic zone con- 
trol effectively using adjacent cell base stations in accordance with the nunr^r of users. 

TO [0003] As the communication speed increases, interference between multipath codes occurs to decrease the com- 
munication quality in an environment such as an urban or indoor space. When the cell area decreases, interference in 
the same channel occurs. 

[0004] One of effective means for removing interference between multipath codes and interference in the same chan- 
nel is to introduce an adaptive antenna effective for removing a long delayed wave and an adaptive equalizer effective 

15 for removing a short delayed wave. In this case, to effectively operate the adaptive antenna and adaptive equalizer, the 
arrival angle and delay time of a signal incoming to the antenna must be measured. Also in dynamic zone control, par- 
ticularly the arrival angle of a signal incoming to the antenna must be measured for beam control of the antenna. 
[0005] Known examples of an algorithm of measuring the arrival angle, delay time, and relative power of a signal 
incoming to the antenna in the multipath environment are a 2D-MUSIC method. FFT-MUSIC method, and 2D-Unitary 

20 ESPRIT method. 

[0006] Evaluation of the indoor propagation environment by the 2D-MUSIC method is reported in "High- Resolution 
Analysis of Indoor Multipath Propagation Sti-ucture", lEICE TRANS. COMMUN.. VOL. E78-B, NO. 1 1 , pp. 1450 - 1457, 
NOVEMBER 1995 (reference 1). 

[0007] The principle of the FFT-MUSIC method is described in "Estimation of Propagation Delay Time and Direction 
25 of Arrival of Indoor Quasi-Millimeter Multipath Waves Using FFT-MUSIC with Triangular Antenna Array", TECHNICAL 
REPORT OF lEICE. A • P95-120. pp. 79 - 84. February 1996 (reference 2). 

[0008] Evaluation of propagation characteristics in the urban space using the 2D-Unitary ESPRIT metiiod is 
described in "High-Resolution 3-D Direction-of-Arrival Determination for Urt>an Mobile Radio". IEEE TRANSACTIONS 
ON ANTENNAS AND PROPAGATION. VOL. 45. NO, 4. APRIL 1997 (references). 

30 [0009] The principle of simultaneously estimating a combination of the arrival angle and delay time of a signal incom- 
ing to the antenna is described in "Paired Estimation of Propagation Delay Times and Directions of Arrival of Multipath 
Waves Using 2D Unitary ESPRir . TECHNICAL REPORT OF lEICE, A • P97-78. JULY 1997 (reference 4). 
[001 0] However, these conventional methods of estimating the arrival angle and delay time of a signal incoming to the 
antenna suffer the following problems. 

35 [0011] First, the 2D-MUS1C method in reference 1 and the 2D-Unitary ESPRIT method in reference 4 can estimate 
only the azimuth plane of the arrival angle and cannot estimate the elevation plane. At present, since a mobile tele- 
phone base station in operation controls the elevation plane directivity of the antenna in accordance with the situation, 
the azimuth and elevation planes of the arrival angle must be simultaneously grasped. 

[0012] Second, the position of a reception antenna element is mechanically moved for measurement. This is also the 
40 same in the 2D-Unitary ESPRIT method in reference 3. The mechanical movement of the reception antenna increases 
the measurement time. Considering mobile communication, the antenna and equalizer must be controlled in real time, 
so these methods are difficult to apply to mobile communication. 

[001 3] Third, since a reception signal is normalized using a transmission signal as a reference, the transmission signal 
must be supplied to the receiver via a cable. Connecting communication stations by the cable essentially contradicts 
45 mobile communication. 

[0014] To the contrary, the FFT-MUSIC method in reference 2 can simultcuieously estimate the arrival angles of azi- 
muth and elevation planes and the delay time without requiring any mechanical movement of tiie reception antenna and 
any connection cable between communication stations. However, this method does not resolve the arrival angle by a 
superresotution technique, resulting in low angle resolution and low precision. In reference 3, although the arrival 
so angles of elevation and azimuth planes are estimated by the 2D-Unitary ESPRIT method, the delay time is estimated 
using a conventional channel counter, resulting in low time resolution and low precision. 

[0015] It is an object of the present invention to provide a radio environment analysis apparatus capable of simulta- 
neously estimating a combination of the arrival angles of elevation and azinrujth planes and the delay time with high res- 
olution and high precision. 

55 [001 6] It is another object of tiie present invention to provide a radio environment measurement apparatus applicable 
to even mobile communication. 

[0017] To achieve the atx)ve objects, according to the present invention, there is provided a radio environment anal- 
ysis apparatus comprising a receiver for receiving a signal modulated by phase shift keying (PSK) using a pseudo noise 
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(PNI).code sequence by an array antenna at a plurality of fixed points on a plane, converting the plurality of received 
signals into intermediate frequencies or demodulated signals, and outputting the resultant signals, and a signal proces- 
sor for processing a signal output from the receiver to analyze an arrival angle, delay time, and relative power of the 
signal incoming to the receiver, the signal processor having a plurality of normalization means for individually normaliz- 

5 ing respective signals from the receiver and outputting the normalized signals as a normalized signal group, first esti- 
mation means for calculating a covariance nr^trix from the signal group nomialized by the normalization means, and 
estimating an arrival angle of each incoming signal using one of a 2D-MUSIC method and 2D-Unitary ESPRIT method, 
arithmetic means for calculating a weight of the array antenna for suppressing incoming signals except for an incoming 
signal having the arrival angle estimated by the first estimation means, and calculating a product of the calculated 

70 weight of the an-ay antenna and the normalized signal group, and second estimation means for calculating a covariance 
matrix from a product signal group output from the arithmetic means, and estimating a delay time and relative power 
using one of the MUSIC method and ESPRIT method. 

Fig. 1 is a block diagram showing a radio environment analysis apparatus according to the first embodiment of the 

15 present invention; 

Fig. 2 is a block diagram showing signal processing operation of a DSP shown in Fig. 1 ; 

Fig. 3 is an explanatory view showing moving average calculation processing by the 2D-Unitary ESPRIT method; 
Fig. 4 is an explanatory view showing a moving average calculation method using the MUSIC method: 
Fig. 5 is a block diagram showing a radio environment analysis apparatus according to the second embodiment of 
20 the present invention; 

Fig. 6 is a block diagram showing signal processing operation of a DSP shown in Fig. 5; and 

Fig. 7 is a block diagram showing a radio environment analysis apparatus according to the third embodiment of the 

present invention. 

25 [0018] The present invention will be described in detail below with reference to the accompanying drawings. 

[001 9] Fig. 1 shows a radio environment analysis apparatus according to the first embodiment of the present inven- 
tion. In Fig. 1. the radio environment analysis apparatus comprises a transmitter 1 for transmitting a signal modulated 
by PSK (Phase Shift Keying) using a PN (Pseudo Noise) code sequence, a receiver 2 for outputting signals received by 
a plurality of antenna elements arrayed on the plane after conversion into intermediate frequencies or demodulated sig- 

30 nals. and a signal processor 3 for simultaneously estimating a combination of the an-ival angles (elevation and azimuth 
angles), delay time, and relative power of a reception signal on the basis of an output from the receiver 2. 
[0020] The transmitter 1 is constituted by a PN code generator 1 1 for generating a PN code, a BPSK (Binary PSK) 
modulator 1 2 for phase-modulating the PN code output from the PN code generator 1 1 . a frequency converter/amplifier 
13 for converting the modulated signal output from the BPSK modulator 12 into a predetermined measurement fre- 

35 quency, and amplifying and outputting it. and a transmission antenna 14 for radiating the signal output from the fre- 
quency converter/amplifier 13 into the space. 

[0021] The receiver 2 is constituted by a planar-array antenna 15 having M x N (M and N are integers equal to or 
larger than 2) antenna elements ISn to 16mn an-ayed on the plane, and a plurality of frequency converters 17ii to 
l^MN for outputting signals received by the antenna elements 16^ and 16^^ after conversion into intermediate fre- 
40 quencies or t>asebands. 

[0022] The signal processor 3 is constituted by a plurality of A/D (Analog-to-Digital) converters (A-D) 1 8i to 1 B^n ^or 
digitally converting the frequency-converted reception signals from the frequency converters 1 7^ ^ to 1 7mn. and a digital 
signal processor (DSP) 19 for simultaneously estimating a combination of the an-ival angles (elevation and azimuth 
angles), delay time, and relative power of a signal incoming to the planar-array antenna 15 on the basis of the plurality 
45 of digital reception signals from the A/D converters 1 S-n to IS^n and a criterion signal. 

[0023] Signal processing operation of the DSP 19 will be explained with reference to Fig. 2. 

[0024] In Fig. 2, a plurality of signals that are received by the antenna elements ISn to 16mn of the planar-array 
antenna 1 5 and digitally converted are transformed by fast Fourier transformers 101^ to IOI^n. Then, the time-series 
signals are converted into signals in the frequency range. To resolve a plurality of correlative signals such as multipath 
so signals for the delay time, a moving average must be calculated in the frequency range. In calculating the moving aver- 
age, the spectrum must be flat in a target frequency range. 

[0025] The group delay devices of the transmitter and receiver must be excluded from processing targets for waiting 

for reception. 

[0026] For this purpose, in an FFT 1 03, a transmission signal is directly input to each channel of the receiver and sub- 
55 jected to fast Fourier transform. In normalization units 102ii to- 102mn. reception signals converted into the frequency 
range are normalized by a calibration (criterion) signal from the FFT 103 that is converted into the frequency range. 
[0027] The arrival angles (elevation and azimuth angles) can be estimated from a reception signal having a frequency 
unique to the planar-array antenna. In a covariance matrix estimation unit 104, a moving-averaged covariance matrix is 
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estimatecJ using the reception signals normalized to the center frequency of the band from the normaliralion units 102,1 
to 102mn- 

[0028] Fig. 3 shows moving average calculation processing by the 2D-Unitary ESPRIT method. 
10029] In Fig. 3. a complex signal received by the planar-aray antenna is divided into (M-m+1) x (N-n+l) , m x n 
5 minor matrices X^. In general, m Is (M+1)/2. and n is (N+1)/2. Each divided reception signal matrix Xpq is vector-trans- 
formed as given by 



10 to calculate a correlation matrix R defined by 



Xpq = vecpCpq) 



[0030] This operation is done for each divided minor matrix Xpq to average a covariance matrix Rpq by 



15 



R = 



1 



(M-m+1) (N-n+l) 



M-m+l N-n+1 



(1) 



(2) 



(3) 



ps1 q°1 
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[0031] The averaged covariance matrix R is unitarily converted by 



(4) 



25 where 



30 



Qn,n=Qm ^Qn 



Qm: mth-order unitary matrix 
Qn: nth-order unitary matrix 
9: Kronecker's operator 

35 [0032] That is, if m(n) is an even number (m(n) « 2k : k is an integer). 

r'k i't 



'2k 



j_r'k i'ki 
^^bik-i'ij 



(4A) 



40 



[0033] If m(n) Is an odd number ( m(n) = 2k+1 ). 



45 



J2 



r'k 0 i'k 



(4B) 



so where 0 is a 0 vector, L is the kth-oider unit matrix, and IIk is defined by 



II, = 
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[0034] To increase resolution and precision, an ensemble average is calculated in a real matrix generation unit 1 05 tsy 

Ry = Re[E[Q"™RQ^]] (5) 

5 

where EQ is the ensemble average. 

The real part of the covariance matrix is extracted. 

[0035] An arrival angle estimation unit 1 06 estimates a combination of the arrival elevation and azimuth angles by tiie 
following procedure. 

10 [0036] First. Ry is decomposed into eigenvalues, and the obtained eigenvalues >^ (i = 1 1 to mn) are arranged in the 
decreasing order to be expressed by 

X^^X^^ ... s » >^L.l = >^mn = (^) 

15 where is noise power. 

From equation (6), the number L of arrival signals can be estimated from the distribution of the eigenvalues. 
[0037] Then, a partial signal space Es to be defined is calculated by 

20 

where ei is an eigenvector corresponding to the eigenvalue X|. 
[0038] Matrices E^ and Ev to be defined are calculated by 



25 Ejj 



where 

30 



=[k::^:]-=[k::::] 



Kv1=K3«>In. K^2=K4®In 
35 K 1 = Re[Q " N-i J nL Kg = lm[Q " ^.^ J nI 

K3 = Re[Q" M-nJaQwl. K4 = lm[Q" M-iJaQwl 

li^: Mth-order unit matrix 
40 l(g: NJth-order unit matrix 

J2: selection matrix for extracting the Mth or Ntii order from the second row of the matrix 

[0039] A matrix Exy to be defined is calculated by 

45 ExY=(KmEs|K^2Es] (9) 

Further, a matrix: 



so 



Exy "Exy = P'^^^'^^^VmEs I K^^Esl (10) 

L(Kji2Es) , 



is calculated and decomposed into eigenvalues. 
55 [0040] A matrix E is calculated by arranging the obtained eigenvectors e. (i = 1 to 2L) by 



E = [eie2...e2j OD 
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The matrix E is divided into four L x L matrices by 



Em 


E|2 




XT 


^21 


^22 



- . . (12) 
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55 



A matrix 4^^ is calculated by 



[0041] Subsequently, this processing using equations (9) to (13) is repetitively executed for to define a matrix 
corresponding to the matrix "V^. "V^ + fVy is decomposed into eigenvalues. Letting the obtained eigenvalue be 
a)^(i) + jco^.(i) (i = 1 to L) . an arrival elevation angle 6j and azimuth angle can be obtained by solving simultaneous 
75 equations: 

to (i)= ^ Ax cos Oicos 4»i (14) 



G)^.(i) = X 0 i sin <t>i 



where Ax and Ay are array element intervals and X is the wavelength of a transmission signal. 
25 [0042] In this way, the arrival azimuth and elevation angles of L arrival signals can be obtained. 
[0043] After that, the delay time and relative power are obtained. 

[0044] An arrival signal designation unit 107 designates arbitrary one of the L detected arrival signals. An adaptive 
weight calculation unit 1 08 stores the designated signal to calculate the weight of the adaptive antenna for suppressing 
other unnecessary signals. 

30 [0045] Letting be the covariance matrix and a{\\f, 9) be a steering vector indicating the direction of the stored signal 
when an unnecessary signal is received by the array antenna, the weight of the adaptive antenna is calculated by 

W=nS^a*(*. e) (15) 

35 An array calculation unit 1 09 calculates the product of the array reception data and weight by 

y = w"x (16) 

[0046] At this time, the product is calculated for a complex reception signal having a plurality of frequencies normal- 
40 ized in the frequency range. If J data are used in the frequency range, array outputs y(l) to y(J) are obtained. Using the 
array output data y(1) to y(J), an estimation moving averaging unit 110 calculates a correlation matrix S by 

S = YY" (17) 

45 where 

Y = [y{i)...y(J)]'' 

[0047] To resolve correlative signals for the delay time, moving averaging shown in Fig. 4 is done for the calculated 
so correlation matrix. Rg. 4 shows a moving average calculation method using the MUSIC method for J = 8. The correla- 
tion matrix S after moving average calculation is decomposed into eigenvalues, and the obtained eigenvalues are 
arranged in the decreasing order to be expressed by 



^ ^2 ^ ... ^ ka » X^j^i = = (18) 



[0048] From equation (18). the number d of arrival waves can be estimated. 

[0049] Then. J-d eigenvectors corresponding to eigenvalues equal to noise power are extracted to calculate a MUSIC 
spectrum P{t) by 
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5 

where 

a(T) = [ej"^V...e'"'^]^ 

70 a>1 ...oxJ: frequencies in use 

[0050] In this evaluation function, a point exhibiting a sharp peak is an estimate of the delay time. The relative power 
is given by 

IS P = (A* A) A*(S - a ^ l)A(A*A) (20) 

where 

A=[a(TO a(Td)] 

so 

I: unit matrix 

[0051 ] Finally, a delay time/relative power estimation unit 1 1 1 calculates an estimate of power from the diagonal term 
of the matrix R Note that processing from the an'ival signal designation unit 1 07 to the delay time/relative power estlma- 
25 tion unit 111 is repetitively executed for all arrival angles of sequentially extracted signals having arrival angles esti- 
mated by the arrival angle estimation unit 106. 

[0052] In the first embodiment, processing of the arrival angle estimation unit 1 06 and processing of the relative power 
estimation unit 111 respectively employ an algorithm based on the 2D-Unitary ESPRIT method and an algorithm based 
on the MUSIC method. Instead, the arrival angle estimation unit 106 and relative power estimation unit 1 1 1 may respec- 
30 lively employ the 2D-MUS1C method and ESPRIT method. 

[0053] Rg. 5 shows a radio environment analysis apparatus according to the second embodiment of the present 
invention. 

[0054] In Fig. 5. a transmitter 1 comprises a frequency sweep signal generator 21 for generating a frequency sweep 
signal as a transmission signal. The transmission frequency from the transmitter 1 changes within a predetermined ^ 

35 band by the frequency sweep signal generator 21 under the control of a signal processor 3. 

[0055] Part of the transmission signal is extracted by a coupler 22 and sent to a receiver 2 via a cable 23. The receiver 
2 additionally comprises a frequency converter 17mni+i for converting the transmission signal input via the cable 23 into 
an intermediate frequency or demodulated signal, amilar to an array antenna reception signal, and outputting the con- 
verted signal to the signal processor. The signal processor 3 comprises a control signal generator 24 controlled by a 

40 DSP 19 in order to change the transmission signal frequency of the transmitter 1 within a predetermined band. 
[0056] Signal processing operation of the DSP 19 shown in Fig. 6 will be explained. 

[0057] Transmission and array antenna reception signals frequency-converted by the receiver 2 are transformed by 
discrete Fourier transformers (DFTs) 1 12ii to 1 12^^ to obtain the amplitudes and phases. The amplitudes and phases 
of the transmission and array antenna reception signals are normalized In normalization units 102^1 to 102^^^. 
45 [0058] The subsequent processing procedure is basically the same as in the first embodiment (see Fig. 2). Note that 
estimating the delay time requires array antenna reception signals having a plurality of frequencies. Hence, array out- 
puts having a plurality of frequencies are calculated In place of the transmission frequency In order to perform process- 
ing of an array calculation unit 109. 

[0059] The second embodiment requires the connection cable 23 between the transmitter 1 and receiver 2, and thus 
so is applicatHe to a static or quasi-static system such as an Indoor radio LAN (Local Area Network). 

[0060] A radio environment analysis apparatus according to the third embodiment of the present invention will be 
described with reference to Fig. 7. In the third embodiment, the present Invention is directly applied to hardware used 
in the above-described references 1 and 4. That is, a receiver 2 comprises a scanner 25 and network analyzer 26, 
whereas a signal processor 3 comprises a personal computer 27. 
55 [0061] As a transmission signal, a frequency sweep signal generated by a frequency sweep signal generator 21 is 
used. The transmission frequency sequentially changes within a predetermined band under the control of the personal 
computer 27 In order to estimate the delay time. Part of the transmission signal is extracted by a coupler 22, input to the 
network analyzer 26 via a cable 23, and used as a reference signal. 
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[0062] A single antenna element 16 is attached to the scanner 25 and controlled by the personal computer 27 so as 
to sequentially move to alignment positions of the anray antenna. The third embodiment obtains an array antenna recep- 
tion signal necessary for estimating the an^ival elevation and azimuth angles by spatial movement of the antenna ele- 
ment 16. 

[0063] The subsequent processing procedure is the same as in the second embodiment (see Fig. 6). Note that since 
the third embodiment adopts the network analyzer 26. the DFTs 112ii to 112^^ and normalization units 102^1 to 
102mn ised to normalize the reception signal can be eliminated. The delay time is estimated using a measurement 
value prepared by calibrating the frequency characteristics of the network analyzer 26. 

[0064] As has been described above, the present invention can estimate not only the anival azimuth angle but also 
the an-ival elevation angle, the delay time, or the propagation parameter of relative power at the same time with high 
resolution and high precision. Further, the present invention can provide a radio environment analysis apparatus appli- 
cable to even mobile communication without requiring any connection cable between communication stations and any 
mechanical movement of the antenna element. 

[0065] The basic processing of the present invention can also be applied to hardware requiring the connection cable 
between communication stations and mechanical movement of the antenna element. 

Claims 

1 . A radio environment analysis apparatus characterized by comprising: 

a receiver (2) for receiving a signal modulated by phase shift keying (PSK) using a pseudo noise (PN) code 
sequence by an array antenna (15, 25) at a plurality of fixed points on a plane, converting tine plurality of 
received signals into intermediate frequencies or demodulated signals, and oulputting the resultant signals; 
and 

a signal processor (3) for processing a signal output from said receiver to analyze an arrival angle, delay time, 
and relative power of tiie signal incoming to said receiver, 
said signal processor having: 

a plurality of normalization means (102ii -1 02mn) ^or individually normalizing respective signals from said 
receiver and outputting the normalized signals as a normalized signal group; 

first estimation means (104 - 106) for calculating a covariance matrix from the signal group normalized by 
said normalization means, and estimating an arrival angle of each incoming signal using one of a 2D- 
MUSIC method and 2D-Llnitary ESPRIT method; 

arithmetic means (108. 109) for calculating a weight of said array antenna for suppressing incoming sig- 
nals except for an incoming signal having the an-ival angle estimated by said first estimation means, and 
calculating a product of tiie calculated weight at said array antenna and the normalized signal group; and 
second estimation means (1 1 0, 1 1 1 ) for calculating a covariance matrix from a product signal group output 
from said arithmetic means, and estimating a delay time and relative power using one of the MUSIC 
mettnod and ESPRIT method. 

2. An apparatus according to daim 1 . wherein said array antenna comprises a plurality of antenna elements (16ii - 
16mn) arrayed on the plane, and 

said normalization means normalizes respective outputs from said plurality of antenna elements. 

3. An apparatus according to claim 1, wherein said array antenna comprises one antenna element (16) moving 
through fixed points on the plane, and 

said normalization means normalizes respective outputs from said antenna element at the fixed points. 

4. An apparatus according to any of daims 1 to 3, wherein said normalization means performs Fourier transform for 
a signal from said receiver and normalizes the signal witii the same spectrum as a transmission signal. 

5. An apparatus according to any of claims 1 to 4. wherein said apparatus further comprises a transmitter (1) for trans- 
mitting a signal PSK-nrodulated using a PN code sequence, and 

said normalization means normalizes tiie signal from said receiver by a frequency sweep signal supplied from said 
transmitter. 

6. An apparatus according to any of claims 1 to 5, wherein said first estimation means estimates a combination of 
arrival elevation and azimuth angles as an arrival angle of an incoming signal, and 
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said arithmetic means stores information of the incoming signal having the arrival elevation and azimuth angles 
estimated by said first estimation means, and calculates a weight of said an^y antenna for suppressing other 
incoming signals in accordance with the storage result. 

5 7. An apparatus according to any of claims 1 to 6, wherein said signal processor stores each incoming signal and cal- 
culates a weight for suppressing other incoming signals to calculate a normalized signal group and product for all 
incoming signals having estimated amval angles, and 

calculates a delay time and relative power from a covariance matrix of the obtained product signal group to esti- 
mate a combination of the arrival angles, delay times, and relative powers of all signals inconrang to said receiver. 

10 

8. A method for analysing a radio environment performing the steps corresponding to any of claims 1 to 7. 
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